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VIBRATION ANALYSIS OF 
CLUSTERED LAUNCH VEHICLES 

SUMMARY 


The Saturn I vehicle represents a complex ensemble of nine beams joined together. 
A first vibration analysis based on an equivalent single beam did not agree satisfactorily 
with the test results of the full scale test vehicle. Comparison of analysis and test re- 
sults rather dictates that the coupling effects of the nine single beam vibrations can not 
be neglected. Hence, the analysis developed in this paper is based on a model consisting 
of nine clustered beams. However, the complexity of the system calls for many simpli- 
fying assumptions. 

The basic idea of the analysis consists in coupling of the single beam vibrations 
such that the connections of the beams will be preserved. The method applied to the un- 
coupled single beams is the well known transfer matrix method. The mentioned coupling 
of the single beam vibrations leads to a linear homogeneous system of compatibility 
equations which must be solved during each transfer. Realization of the boundary con- 
ditions leads to a frequency curve indicating the coupled frequencies of the system. 

Finally the coupled mode shapes can be determined. 

The analysis is programmed on IBM-709-4. Because it consists of simple 
matrix operations the computer time is short. 

The results of an analysis of the Saturn I vehicle compare well with test results. 


SECTION I. INTRODUCTION 


The Saturn vehicle as shown by Figure i consists of payload, S-IV stage and S-I 
stage. In the S-I stage, eight clustered tanks (four fuel and four lox) are attached at 
each end to the center barrel (lox tank) . Then, the Saturn vehicle represents a highly 
complex ensemble of nine beams joined together. 

The results of a first vibration analysis based on an equivalent single beam did 
not agree satisfactorily with the test results of the full scale test vehicle (Ref. i) . The 
comparison of analysis and test results dictates that coupling can not be neglected; indeed, 
the vibrations of the nine single beams must be coupled. 



Subsequently, several proposals concerning the vibration analysis of clustered 
booster vehicles have been made. Some of them are suggested by References 2, 3, and 4. 
The analysis shown herein is that done by the Dynamics and Loads Branch (Structures 
Division of the Propulsion and Vehicle Engineering Laboratory) George C. Marshall 
Space Flight Center. 

This analysis is of a model outlined in Section III of this report. Complexity of 
the Saturn vehicle calls for many simplifying assumptions. But, remember the funda- 
mental premise: the basic configuration of the model is a central beam with eight 
smaller beams clustered about it. From this premise, simplifications can be made 
concerning the mode of vibration, and the support conditions of the outer beams and 
single beam models. Assumptions about the mode of vibration are: 

a. The center beam is restricted to plane vibrations in a symmetry plane of the 
vehicle through its long axis. 

b. The outer beam vibrations occur symmetrically to this plane. 

Nevertheless, each outer beam may have two planes of vibrations. The outer beam 
supports are assumed equivalent to either simple supports, or to simple supports having 
additional bending restraints. The single beams, finally, are modeled by lumped mass 
systems based on Timoshenko's beam theory. 

The essential part of the analysis (Section IV) deals with the coupling of the 
single beam vibrations. The beams do not vibrate independently — they are connected. 
Coupling means adjustment of outer beam and center beam vibrations such that the 
connections of the beams will be preserved. In other words, certain conditions of com- 
patibility must be satisfied. These circumstances can be expressed mathematically by 
a linear inhomogeneous system of equations to be solved. 

The. fundamental method applied in this analysis is the transfer matrix method. 

This method is well known (Ref. 5) but, for completeness a short explanation is given in 
the appendix. The transfer matrix method is based on the assumption that the beam 
masses and mass moments of inertia are lumped at stations of the beam axes. The 
state at each station is determined by a state vector which contains deflection, rotation 
angle, moment and shear at this station as components. The state vectors of successive 
stations are linked by transfer matrices. Such a transfer matrix is the product of an 
inertia and a stiffness matrix. In this way the transfer from one station to any other 
station of the beam is determined by a matrix which is composed of inertia and stiffness 
matrices. This applies especially to the end stations of the beam. These beam transfer 
matrices of the center and the outer beams are defined in Table 2. From the elements of 
these matrices the above mentioned compatibility equations can be formed and finally 
solved. Then the transfer of the complete system can be performed (Table 3). Now, 
realization of the boundary conditions leads to the eigenvalue equation. Evaluation of this 
equation and determination of the mode shapes form the final steps of the analysis. 
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The whole procedure outlined above will be done by application of simple matrix 
algebra. Formation, multiplication, inversion of low order matrices are the only opera- 
tions performed. The matrices used and their dimensions are compiled in Table 4. The 
analysis is programmed on IBM-709-4. Because it consists of the mentioned simple 
matrix operations, the computer time is very short. To obtain one mode shape and fre- 
quency, approximately three minutes are necessary. 

The results of an analysis of the Saturn I, SA-5 vehicle are presented and compare 
well with the test results. The mode shape characteristics of these complex shapes in- 
dicate the proper coupling model is used in the analysis. 


SECTION II. MODEL 

Before beginning discussion of the model used, a short description of the vehicle 
will be given (Fig. 1) . The vehicle consists of payload, second stage and first stage. In 
the first stage, eight tanks are clustered symmetrically about a center tank. To connect 
outer tanks with the center tank, a spider beam at the upper end and the outrigger assembly 
on the lower end are used (Fig. 2) . In addition, spider beam and outrigger assembly 
provide interconnection between stages and supports of the engines, respectively. The 
eight outer tanks consist of four lox tanks and four fuel tanks . The center tank is a 
lox tank. 

The basic configuration of the model is, therefore, a center beam with eight 
smaller beams clustered symmetrically about it. Spider and outrigger structure provide 
supports of the outer beams. 

The preparation of a vibration analysis of this system represents a difficult pro- 
blem. The elastic, vibrating structure consists of a three-dimensional ensemble of non- 
uniform beams connected by the complex spider and outrigger structure. Obviously, the 
analysis must be based on a simplified model. On the other hand, analytic results based 
on too-simple models do not agree satisfactorily with the test results of the full scale 
test vehicle. Studies made, using different models, indicate that realistic assumptions 
concerning the outer beam supports are significant. Thus, to build up a mathematical 
model, two basic principles should be considered: 

a. For analytic reasons a simplified model must be used. 

b. To represent satisfactorily the dynamic behavior of the actual vehicle, the 
outer beam supports of the spider and tail sections must be modeled as closely as possible. 

Some of the simplifying assumptions are usual and well known: 

a. The nine beams of the vehicle model will be considered as Timoshenko beams. 
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b. The beam masses will be lumped at stations of the beam axes such that the 
stiffness between stations may be assumed as constant. 

c. The liquid within the tanks will be assumed as solid — having mass, but no 
mass moment of inertia. The masses will be lumped like the other solid masses on the 
beam axes. 

Another assumption deals with the vibration itself. Obviously a three-dimensional 
cluster, as in the Saturn vehicle, may vibrate in different modes depending on the nature 
of the excitation. For simplification, the analysis proposed is restricted to vibrations 
excited in a plane of symmetry. Cases (a) and (b) of Figure 3 show the positions of 
such planes within the cluster. 1 . Then, the center beam vibration occurs in the symmetry 
plane. The outer beam vibrates symmetrically to this plane. Although the center beam 
vibration is in one plane, the outer beams must be expected to vibrate in two planes. 

The analysis under discussion is based on case (a) of Figure 3. Yet, the analysis 
can be extended to case (b) . On the .other hand, test results of the full scale test vehicle 
do not indicate considerable deviations of mode shapes and frequencies of the cases (a) , 
(b) , and (c) shown by Figure 3. 

The two planes of vibration of the outer beams can be assumed as radial and 
tangential planes (Fig. 4a). To distinguish among the different vibration planes, these 
planes will be numbered ( Fig. 4a) . For reasons of symmetry, restriction to one side 
of the vehicle is possible. The planes of radial vibrations are: 

(1) (2) (4) (6) (8) 

while tangential vibrations occur in the planes: 

(3) (5) (7) 

To measure deflections and rotation angles, coordinate systems must be assigned 
to each vibration plane (Fig. 4b) . Hence the outer beam vibration planes may be 
characterized by the angles a. = i, 2, ... (8) between the y-axes of the center beam and 
outer beam vibration planes. Since only the cosine of these angles is important, the 
angles may be measured in arbitrary directions (Fig. 4a, Table i). 


i. It must be mentioned that these planes do not represent perfect symmetry planes 
since the arrangement of the outriggers is not symmetric to these planes ( Figs. 4 and 5) . 
However, the differences between both types of outriggers are so small — especially with 
regard to stiffness — that the outriggers can be considered as equal. 
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Now, the most important assumptions have to do with the outer beam supports 
located on spider arms and outriggers (Fig. 2). Spider arms and outriggers will be 
regarded as massless beams rigidly connected with the center beam at two branching 
points: Branching point i in the spider section, branching point 2 in the tail section. 

The flexibility of spider arms, and outriggers can be determined by the influence co*- 
efficients <5^, <5^ respectively of the support points. The 6's are the deflections per 

unit load in direction of the tank axes. In view of the outer beam connections shown by 
Figure 2, the following assumptions can be made. 

Spider section 

In the radial planes of vibration (1), (2), (4), (6), (8), and in the tangential vibration 
plane (5) the outer beams behave like simple supported beams. 


Hence, the boundary conditions are 


yi 
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1 

* 

M* 

- cpi cos a.) ; i = 3, 7 


( 2 ) 


where <P\ are the rotation angles of the i^ outer beam and the center beam at 

branching point i respectively. 
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Figure 5 shows a schematic sketch of the support conditions of the tangential 
vibration planes. For convenience, the outer beam axes (in reality, situated on a cylinder 
around the center beam) are rolled up in the drawing plane. One concludes that the 
spring stiffness 

d 2 
o 


where d^ represents the outer tank diameter. 

Since bending of the center beam causes rotation of the spider beam plane about 
(3) ( 7 ) 

the angle (Pi, the angles (p j and <Pi are not the effective angles acting on the 
springs. Rather the effective angles are given by the parenthetical expression of the 
last two equations of the system 2. 

Tail Section 


In the radial planes of vibration (1), (2), (4), (6), (8) ; the outer beams can be con- 
sidered again as simple supported beams. Hence, the boundary conditions of these 
planes may be written: 


(i) 

y 2 = y 2 cos a. 

M«)-0 ‘ < 3 > 

where y 2 ^ , M 2 ^ represent deflection and moment of the i^ outer beam vibration at 
the tail support and y 2 is the center beam deflection of branching point 2. 

The beams vibrating in the tangential planes (3), (5), (7), may be regarded as 
supported beams having additional bending restraints. So, with regard to the deflections 
it follows 


y2 


(i) 


y 2 cos a. ; i = 3, 5, 7 


(4) 


To obtain the remaining boundary conditions of the vibration planes (3) , (5) , (7) , 
Figure 5 may be considered. It follows easily that 
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(5) 


<?2 


( 3 ) 


<P 2 
<P 2 


(5) 

( 7 ) 


<p 2 cos a 3 
<P 2 cos a 5 
cp 2 cos Q!7 


V /d o 

" v m}^ d o 

(^ni ” 77 iv 


O 


l 


where is the rotation angle of the i outer beam vibration at the tail support and 

cp 2 represents the vibration angle of the center beam at branching point i. 17. is the 

relative deflection of the support point i (i = I, II, III, IV) . The terms on the left side of 
equations (5) represent the relative outer beam rotation angles at the support points 
which affect the existing springs. 


Obviously, bending of the center beam can cause axial forces in the outer tanks. 
Still, from Figure 2, one realizes that axial forces may attack the outer lox tanks only. 
The outer fuel tanks may not resist axial forces because of their special supports in the 
spider section. Under consideration of the axial lox tank forces S 3 S 5 it follows from 
Figure 5 that 



Now, from Figure 5, 


it may be concluded 


(<P 2 - <Pi) e 



+ 6 l + 


5 

o 

2 


) 


where e is defined by 


Figure 5 and <5 

J-l 


is the lox tank influence coefficient. 


Setting 


( 6 ) 



= 6 


7 



one obtains 


S3 - g (V2 ~ •Pi) 
and similar 

S 7 = - ~ (<p2 ~<Pi) 

From (5), (6) and (7) it follows 


( 7 ) 


CO 

CM 

- cp 2 COS £*3 = 

n/r ( 3 ) 

C 33 m 2 

+ 

C 35 m 2 (5) 



- <p2 COS £* 5 = C L (<p 2 

- <p x ) + C 53 M 2 (3) + C 55 M 2 (5) 

+ 

(7) 

C 57 m 2 

> (8) 

ft™ 

- cp 2 cos Oil ~ 

C 75 m 2 (5) 

-1- 

A(r ( 7 ) 

c 7 7 M 2 y 



where 


ed 

o 

C L d <5 


(9) 


26 

o 

c 33 _ c 53 “ c 77 “ j 2 


c 35 “ c 53 ~ " ' 


d 2 
o 


> 


( 10 ) 


c 37 _ c 75 “ 


d 2 

O 


Equations 3, 4, 8, 9, and 10 represent the complete boundary condition set of the 
outer tanks at the tail section. 


At this point, some remarks will be made about coupling between the outer beam 
and the center beam vibrations. Obviously, these vibrations can not occur independently 
The center beam distortions of branching point 1 and 2 influence the boundary conditions 
of the outer beams, while the outer beam support forces, and moments, contribute to 
shear and moment of the center beam at branching points 1 and 2. 
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V 


Coupling exists also between the outer tank tangential vibrations. No direct 
coupling exists between the vibrations of radial planes (1), (2), (6), (8). Hence, the 
fuel tank modes of the planes (i) , (8) and the lox tank modes of the planes (2) and (6) 
are equal to each other. (See Figs. 4 and 6) 


In the above boundary conditions (1), (2), (3), (4), (8) , the terms containing 
cos a. (i = 1 , 2, . . . 3) and c imply the coupling of the outer beam vibrations with the 

1 Ju 

center beam vibration. The terms with the coefficient c„, i, j = 3, 5, 7 imply the 
coupling of the tangential plane vibrations (3) , (5) , (7), with each other. 


Because 


cos aq = cos 90° = 0 

the vibration of the radial plane (4) is uncoupled and so, does not contribute to the center 
beam vibration. Conversely, a vibration in plane (4) can not be excited since plane (4) 
is perpendicular to the assumed plane of excitation. For that reason, the radial vibration 
of plane (4) could have been omitted from this analysis without any loss of generality. 


SECTION HI. NUMERICAL RESULTS AND COMPARISON WITH TEST RESULTS 


It is desirable to compare the results of the analysis with some dynamic test 
results in order to substantiate the analysis procedures and the mathematical model to 
the extent that sufficient structural detail was considered. The comparison can best be 
evaluated if the results of the analysis are considered first. The analysis, as described 
in this report, includes some interesting coupling of the outer tanks to the main beam. 

The dynamic test results which are available, also need some explanation since the test 
was a determination of transfer functions rather than a model survey. 

Some typical mode shapes from the analysis are shown in Figures 7 to 16. (See 
also Fig. 6) Each figure represents a single mode of the total vehicle and can be identi- 
fied, with some reservations, as either the center line of the vehicle bending in a first, 
second, or higher beam mode shape; or the bending of one of the outer tanks representing 
the major bending. The modes identified by vehicle center line are called "bending 
modes" and those identified by outer tanks are called "cluster modes. " This identifica- 
tion can be rather useful in the evaluation of the analysis but can also offer some confusion 
in the proper presentation of results. The identification is maintained because of its 
usefulness in analysis and because of the similar identification of the dynamic test 
results. 
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The complication of a shape identification method can be seen if the mode shapes 
are studied as the fill level is lowered in the first stage of the vehicle to represent change 
of propellant with flight time. The simple and rather logical way to show the frequency 
trend throughout flight is to join the lowest frequency with a smooth curve, the next 
lowest, and so forth for the available data. As the total weight is reduced, the expecta- 
tion would be for an increase in the frequency of the system. The frequency trend plotted 
in Figure 17 shows that this is true when these next lowest frequencies are connected. 
When sufficient fill conditions are used, the frequency trend lines will represent the 
frequency trends of the system. 

Now the same data points can be joined by lines which represent the same mode 
shape trend from one time point to another. The trends of the "bending" modes and 
"cluster" modes are shown in Figure 18. The points on these curves between data points 
do not correspond to vehicle frequencies, with the exception of the line segments that 
happen to be identical to Figure 17. When this identification is used, the resulting curves 
can have a decrease in frequency as the fill level is lowered. This decrease does not 
show in Figure 18, but is not unusual. The curves also intersect and cross each other in 
a manner which can be described as the cluster modes crossing the center line bending 
modes. This crossing of the trends of Figure 18 brings up the question of crossing of 
the trends of Figure 17. These trends can cross, or come together, since cases can 
happen, and do happen where two modes of the system cannot be separated and a single 
frequency exists for two different mode shapes. 

Now the test results will be considered. The dynamic testing is conducted by 
hanging the vehicle on a coil spring-cable suspension system in a special structural steel 
test tower. The cables are attached to the base of the booster stage and to the test tower 
at about the level of the top of the booster. A small angle is allowed between the cable 
from the vertical, and the springs between the cable and tower are as soft as practical. 
This system provides little lateral restraint and only a small end moment, which is 
necessary to maintain the vertical stability of the vehicle. 

With the vehicle suspended in the tower, a single shaker is used to conduct the 
test program. The shaker is attached at the engine gimbal station of the first stage and 
transfer functions are determined for the one shaker location. The amplitude at the 
frequencies with maximum response is plotted and decay damping at these frequencies is 
determined. From these tests, the frequency trend curve and the "response mode 
shapes" will be used for comparison. The comparisons shown are as complete as 
practical and a few comments are in order. 

The comparison of frequency trends is shown in Figure i9. The first bending 
mode frequencies compare favorably. The largest deviation is in the second mode and 
can be partly accounted for by the absence in the analysis of an engine mode representation. 
This has been added to later analysis and provides a better comparison. The test did 
not find an engine mode at all test points and is an example of the difficulty in obtaining 
all modes on the dynamic test vehicle. 
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The frequency comparison without mode shape comparison is not sufficient for 
such a system, therefore, the mode shapes are presented for the flight time of 35 seconds. 
The calculated mode shapes show good comparison to center line bending of the test modes. 
Some test mode shapes are compared with two calculated modes, since the test is a 
"response mode shape." The outer tanks do not compare as well in amplitude as the 
center line, since the tank response is generally lower. The comparison to more than 
one calculated mode shape and the difference in outer tank relative amplitude can partly 
be attributed to the type of test conducted. 

The results from the analysis give good agreement with test results. It can be 
concluded that the structural connections and general structural parameters have been 
properly used since the mode shapes compare favorably with test results. The structural 
spring constants for tank coupling have been improved since this analysis and theoretical 
response analysis using the calculated mode shapes indicate the lower amplitude response 
of outer tanks is to be expected and the response of adjacent modes can be significant. 


SECTION IV. ANALYSIS 


The transfer matrix method as outlined in the appendix will now be applied to the 
developed model. The notations of the state vectors and transfer matrices are given in 
Tables 2 and 3. The simple transfers of the single beams (center beam, outer beams) 
may be seen from Table 2 while the transfer of the complete model is given by Table 3. 
The inertia and stiffness matrices T and S will be determined from the given vehicle 

( F) ( L) 

data. (See appendix. ) For this reason, the matrices Lj, L 2 , L 3 , L 2 , L 2 

which are composed of inertia and stiffness matrices (see appendix) can be considered 

as known. 

From Table 3 of the transfer O— 3 of the model may be concluded as 


y 3 *-3 ^2 *-2 *“1 ^1 y 0 - 1-3 Ljj'Lj y 0 - L y„ 

l* = l.;l 2 l; 

- L I | 

S “ L 3 L 2 L 1 


where 


Li = 


(i) 

l ( 1 ) 
*12 

i (i) 

*13 

i d) 

n 

*14 

(i) 

( (i) 

i (i) 
* 23 

# ( 1 ) 
*24 

21 

*22 

( 1 ) 

p (!) 

,<i) 

*33 

I ( 1 ) 

31 


*34 

(i) 

M) 

*42 

, ( 1 ) 
*43 

t ( 1 ) 

*44 

41 


i = i, 2, 3 


(ii) 


( 12 ) 
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L' t represents the transfer i — 1' (over branching point i) and L' 2 represents the transfer 
2' — 2 (over branching point 2) . 

To perform the transfer 0 — 3 the matrices L'* and L 2 * must be known. The 
determination of these matrices represents the main problem of the analysis and will be 
done in two steps. First, will be determined and then in a direct way L 2 * . L' 2 will 
not be determined explicitly. 

The sketch on the top of Table 3 shows the vehicle model represented by the 
center beam axis. Because the outer beams are removed, external forces and moments 
must be applied to ensure equilibrium. 

It follows, then, from this sketch: 


y'i = y i + 


y2 - y 2 + 


0 

o 

AMi 

AQi 


0 

0 

AM 2 

aq 2 


(13) 


(14) 


AQ^, i = 1, 2 is the resultant of all outer beam support forces. One may, 


therefore, write: 


AQ. = (-1) 1 [Q. U) cos cos o; 8 + 2^ Q.^ cos a.] ; i - 1, 2 

1 1 1 j=2 1 J 


,(j) 


where 


; i = 1, 2; j = 1, 2. .. 8 are the shears at the outer beam supports. 


Using vector and matrix notation, one may write: 
= (-1) 1 V 4 q.; i = 1, 2 


0 

0 

0 

AQ. 


(15) 


12 



where 


0 

0 

0 

COS C*! 


0 

0 

0 

2 cos a 2 


0 0 

0 0 

0 0 

. 2 COS Q! 7 COS CHg 


(16) 


and 


q, = 


Q i 

Q; 


( 1 ) 

( 2 ) 


S (8) 


i = 1, 2 


(17) 


AM^ is caused not only by the outer beam support reactions, there is also a con- 
tribution of the lox tank axial forces. AM^ can, therefore, be written: 


AM. = AM! + AM'.' ; i = 1, 2 

1X1 


(18) 


where 


AM! 1 represents the moment contribution of the outer lox tank axial forces. 
AM^ can be expressed similar to equation 15 


0 

0 

AM! 
0 1 


= (-i) v 3 m. 


(19) 


i = 1, 2 


where 


m . = 


M 


(i) 


to 


m; 

1 


m (8 > 

1 


(20) 


i = 1, 2 


and V 3 results from the matrix 16 by changing third and fourth row. 
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The second term of 18 may be obtained from Figure 5 and equation 7 as 
AM." = 2(-l) i + 1 2eS 3 = (-l) i+ 1 (<P 2 - <P\) 


or 


0 

0 

AM'.' 


= (-i) 1+ 1 c u* 4 (y'i - y[) 


i = i, 2 


where 


c = 


4e 


and 


U*, = 


0 0 0 0 
0 0 0 0 
0 i 0 0 

0 0 0 0 


In equation 21 y 2 may be replaced by 

y 2 - l 2 y i 

(Table 3) then it follows 


0 

0 

AM'. 

l 

0 


(-1) 1+1 c A yj 


i = 1, 2 


> 


f 


A = U 4 (L 2 - U 4 ) 

where U 4 is the four x four unit matrix. 

From equations 13, 15, 19, 24 it follows: 

y\ = Yi - v 4 - v 3 m i +c A y\ 


( 21 ) 


(■ 22 ) 


(23) 


(24) 


(25) 
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Now to replace m, the boundary conditions ( 1 ) and (2) with regard to 
j = i, 2 8, may be introduced. These conditions, written in vector notation are: 


= - K f ! + K y t 


( 26 ) 


where 


K = 


0 0 0 
0 0 0 


0 K o 0 0 0 


0 
0 
0 

0 0 0 0 0 0 0 


0 0 0 0 
0 0 0 0 
0 
0 


00000000 

00000000 

000000/C0 

00000000 


(27) 


f. = 

X 


<p 


(1) 


5 2) 


*i 




i = 1, 2 


(28) 


V * 2 = 


COS a* 

cos a 2 
cos a 3 

COS Q! 4 

cos a 5 
cos a e 
cos a 7 

COS Qfg 


0 0 
0 0 
0 0 
0 0 


0 

0 


0 0 
0 0 


(29) 


Insertion of equation 26 into equation 25 — and little algebra — yields 
(U 4 - C A) y ; = (U 4 - V 3 K V* ) yi + V 3 K fx - v 4 q t 


(30) 


Setting 


(U 4 - c A) 


= U 4 +• 


A =R 


i l< 2 > 

• - C 1 23 


(31) 
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and 


R(U 4 -V 3 KV*)= B t ] 

R V 3 K t = - B 2 > 

R v 4 = b 3 J 

equation 30 changes over in 

yi = BiVi - B 2 f 1 - B 3 q t 

The relation 3i may be proved easily by verifying 


(32) 


(33) 


(U - c A) (U +— ! £ A) =U 

' - C I 23 

From equations 14, 15, 19, 23, 24 one obtains 

y 2 = (L 2 - c A) yj +v 3 m 2 + V 4 q 2 (34) 


Equation 34 will be needed later on. 

In equation 33 the sixteen components of the vectors and q 4 represent the un- 
knowns of the problem in hand. In the following it will be shown that f 4 and q 4 can be 
determined by the outer beam boundary conditions of the tail section. From these con- 
ditions it will be concluded: 

;* :V* } 

'll - Q y 4 J 

where F and Q are eight x four matrices. 


Then, combining equations 33 and 35, one obtains 

y ;=(Bi-B 2 F-B 3 Q) Yl (36) 

and hence (see Table 3) 

i-l =B l ‘ B 2 F * B 3 Q ( 37) 
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To determine the matrices F and Q of equations 35 the transfers of the outer beams 
must be considered. In doing so, it is convenient to use the following notations: 


^ il 


y* 


(i) 


U 


y« 




’ ^i2 - ^i' ^ i3 m i' ^ i4 ^i 


i= 1, 2 


where q^, rrr, f. are defined by equations 17, 20, and 28, respectively. 
Then the outer beam transfers (Table 2) may be written: 


(38) 


Z 1 L kj d t y k = l, 2, 3, 4 


(39) 


where 


L kj 


,(F) 

kj 


,(L) 

kj 


0 

0 

0 

0 

0 

0 


,<L) 

'kj 


0 

0 

0 

0 

0 


0 

0 

0 


,(F) 

'kj 


0 

0 

0 

0 


0 

0 

0 

0 

4 

0 

0 

0 


i ( F) 

kj 


0 

0 

0 

0 

0 


, (L) 

'kj 


0 

0 


0 

0 

0 

0 

0 

0 


,(L) 

kj 


0 

0 

0 

0 

0 

0 

0 


i ( F) 

kj 


k,j = 1, 2, 3, 4 


(40) 


Application of the notation 38 to the boundary conditions: equations 1, 2, and 26, 
respectively, yields 
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d il = v l' Y1 


d 13 = m i -•Kf 1 +KV 2 ! 'y 1 


( 41 ) 


where K and V 2 * are given by equations 27 and 29. Vi* follows from the matrix 29 by 
changing the first two columns. 


where 


Equations 39, 40, 4i result in 
d 2i- = D iYi + E i f i + L-i4 

Di =Lit V^ + Lis KV? 

E i ~ E i2 ‘ L i3 K 






i = i, 2, 3, 4 


(42) 


J 


As already mentioned, the boundary conditions 3, 4,8, will now serve for deter- 
mination of fi and qj. Using the vector notation 38 these conditions may be written: 


d 2 , = v? y ; 


u8*< d 22-<y;> = c L ( y ;. y ;) + c M d 
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} 


(43) 


where 


U f 


0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

o" 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

; C L = 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 


0 

C L 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 
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i 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

C 33 

0 

C 35 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

C 53 

0 

C 55 

0 

C 57 

0 

0 

0 

0 

0 

0 

i 

0 

0 

0 

0 

0 

0 

C 75 

0 

C77 

0 

0 

0 

0 

0 

0 

0- 

0 

1 


The elements C T ; C..;i,j=3, 5, 7 are given by equations 9 and 10. 
1j ij 

Insertion of equation 23 into equation 43 and collection of terms yields 


d 2 l = Vf L 2 y{ 

U 8 ^22 


C M ^23 CYi 


(44) 


where 


C = (U e V 2* + Cl) L 2 - C L 




Now, replacement of y] , d 2 i , d 22 > <J 2 3 by equations 33 and 42 changes the first two 
equations 44 over in 


N 11 +N 12 f l " R 1 Vl 1 

N 21 9 1 + N 2 2 ^2 _ ^2 y 2 J 

where 

Nu = L 14 + V! ! L 2 B 3 
N i 2 =E 1 +VfL 2 B 2 
Rr-Di + vt L 2 B ! 

N 21 = U 8 L 2 4 +CB 3 -C ]y[L 34 
N 22 = Us E 2 +C B 2 * C m E 3 
R 2 = * U 8 D 2 +CB! + CjyjD 3 


(45) 


(46) 
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The matrices N i, j = 1, 2 are eight x eight matrices, while R t R 2 are eight 
x four matrices. Solution of this system is given by equations 35 where 


F = r 1 [R t - N u Ng^R,] 

Q = N 2i [- N 22 I 1 R i + N 22 ^ ^ 

I - N 12 - Nn Ngj 1 N 22 

Proof by insertion. 


(47) 


Obviously the system of equation 45 is solvable only if the matrices I and N 2 i are 
nonsingular matrices. 

Now from equations 37 and 47 it follows that Lj is known. 

To determine L 2 * given by the second equation 11 and Table 3, one may start 
from equation 34. Considering equations 33, 35, and 42 it follows from equation 34: 

y 2 = (J + g F + H Q) yi = L 2 ' yi (48) 

where J = V 3 D 3 +V 4 D 4 +(|_ 2 -cA)Bi 

9‘ =V 3 E 3 +V 4 E 4-(L 2 ^ A ) B 2 l 

f (49) 

H = V 3 L 34 + V 4 L 44 - (L 2 - c A) B 3 J 

Because of 47, 48, and 49 L 2 * is known and hence 

L = L , L »*L 1 

also (See eq. 13) . In this way, the frequency equation of the system can be solved as 
outlined in the appendix. The solutions are the eigenvalues 

*1 = «i ; A 2 = 0)| . . . 

Now the start vectors y can be determined. 

7 o 

y is the start vector of the center beam part between top 0 and branching point 1 . 
The start vector of the center beam between branching points 1 and 2 is given by 

L ! L i 7o 
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The start vectors, vectors of the outer beam vibrations, are 


where 


yi 

<p i 

Mi 


(i) 

(i) 

(i) 

<i> 


_Qi i i = 
The components 


y 1 = L i y ( 


i, 2, ... 8 

of these start vectors follow from equations 35, 4i, and 47, 
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FIGURE 2. SATURN I BLOCK II CONNECTION DETAILS 
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FIGURE 4. COORDINATE SYSTEMS 
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FIGURE 6. SIGN CONVENTION FOR THE VEHICLE MODE SHAPE PLOTS 

The sign convention for the vehicle mode shape is positive as indicated by the 
arrows in the Figure above. The deflection curves of the tanks are plotted considering 
the following: 

1. The deflection curves are plotted in the planes as indicated above and are 
represented in the vehicle plane of symmetry. 

2. The reference axis for all deflections is the undisturbed main vehicle axis. 

3. The eight planes are plotted with symmetry considered. Planes 1, 2, 3, 5, 
and 7 are plotted with their sign convention. Planes 6 and 8 have the same deflections 
as planes 2 and i, respectively, with opposite signs. Plane 4 has no deflection. 

Notice: This sign convention is opposite to that used in the analysis. 
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FIGURE 12. SATURN I (SAD-5) VEHICLE, RELATIVE AMPLITUDE VERSUS VEHICLE STATION 
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FIGURE 13. SATURN I (SAD-5) VEHICLE, RELATIVE AMPLITUDE VERSUS VEHICLE STATION 
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FIGURE 14. SATURN I (SAD-5) VEHICLE, RELATIVE AMPLITUDE VERSUS VEHICLE STATION 
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FIGURE 15. SATURN I (SAD-5) VEHICLE, RELATIVE AMPLITUDE VERSUS VEHICLE STATION 
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FIGURE 19. SATURN I (SAD-5) FREQUENCY COMPARISON OF REPRESENTATIVE 
MODE SHAPES 
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Table 3. Transfer of the Center Beam 




AQ 1 


Center Beam Axis 




AM 2 

“3 


Total Launch Vehicle 
Model 


Top 

0 


1 

l 

1' 

1 

2' 

1 

2 

1 

End 

3 

State 

Vectors 

yo = 

| 


1 

Yi = 

Yi 

<Pi 

Mi 

Qi_ 


yi = 

yi 

<?i 

Mi 

_Ql_ 


t 

y 2 = 

y2 

(p 2 
Mi 

q 2 

1 

YZ = 

y2 

Wz 

M, 

q 2 

L. — 


Yz = 

y% 

<?3 

Mg 

Q 3 


Transfer 

Matrices 


Li 


Lj 




l 2 




Transfer 
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Yi =Li y 0 


! 1 T 
yi = L iyi 


• 1 1 

Yz = L 2 yi 


l ' ’ 

y 2 = L 2 y 2 

□ 

ya - L 3 y 2 





























Table 4. Used Matrices and Their Dimensions 


Dimension 

Matrices 

four x four 

L., i = i, 2, 3 ; Lj, L 2 , L 2 * , L, L 2 ^ , L 2 ^ 
Bi, J j To i - i, 2 , . . . n, L,Ej, 

four x eight 

V4, V3, &2> B3, G, H 

eight x four 

V 2 * > F, Q, Vi*, D r i= i, 2, 3, 4, C L , R r 


K, L„ ,i,j — 1 , 2, 3, 4 ,Ej, i — i, 2, 3 , 4 , U| 

i,j = 1, 2, 1 

eight x eight 


Ur > U4> A, R, 


i= i, 2 


’ C M’ N ij’ 
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APPENDIX 


TRANSFER MATRIX METHOD 
APPLIED TO A LUMPED MASS BEAM 


A spring-supported lumped mass beam is shown by Figure 20a. It consists of masses 
m. and mass moments of inertia P. concentrated at Stations i on the beam axis (i = i, 2, 

3 n) . The stiffness between successive stations is constant. 

If it is assumed that this system is free-vibrating, following a sinusoidal law, 
shear and moment will obey this law also. The equilibrium conditions of the i^ mass 
(Fig. 20b) written in the amplitudes of the mentioned quantities are 


M ,* = X p. q> . + M . 

l ill 


> (Ai) 

V =- Xm i>'i J 

where y., are deflection and rotation angle at station i, m., q.. Q* represent 

th 

moment and shear at left and right side, respectively, of the i mass point. X is 
the unknown parameter, of the problem, 


X = co 2 = 4 7T 2 f 2 (A2) 

where a> is the circular frequency and f is the natural frequency of the vibrating system. 
The state in a point of the lumped mass system is determined by the state vector, 
y 

<p 

y " M 

Q 
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From equation Ai it follows 


y -f 


1 0 0 0 


y. 

1 




i 

(P ^ ^ 

= 

0 i 0 0 



M.* 

i 


0 AP. i 0 

l 


M. 

1 

Qj* 

k. - 


1 

1 

r 

o 

o 

1 


Qi 


(A3) 


or y.* = T y ( A4) 

1 1 i 

where the "inertia" matrix Tj is defined by equation ( A3) . 

A similar relation exists between the state vectors y.* and y (Fig 20c) . 


y i+i 



(A5) 


where S is the stiffness matrix. S ^ can be formed by solving the linear system of 
differential equations: 


GA si (y' - 9) = Q \ 
- El 9' = M J 


(A6) 


with the boundary conditions 


x 


y-7ii 


9 = 9j 


(A7) 


where (Fig. 20c) 


M = M * + (x-x) Q .* 
X 11 

Q = Q * 

l 


(AS) 
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The solution of equations A6 (x = x. + is given by 


y i+ 1 " y i + l * * 


l } /£ , 3 

1 M - 1 


i i 2EI . i 

l 


<p 

i.+ i 


£ . 1.2 
1 M * - * 


6E I'. GA 
‘ i si. 


9 i _ El. V ” 2eT V 

l l 


> 


(A9) 


where £ . = x. , . - x. 

i l+l i 


From equations A5, A7, A8, and A9 it follows 


S i = 


1 £. -!?/ 2EI. 

l li 


0 0 
0 0 


0 


(!f/6EI. - l./GA .) 
i l l si 7 


0 1 -l./EI. - 1 ? /2EI 

ii l i 

1 


i. 


1 


Repeated application of equations A4 and A5 results in 


y = S T y 

2 l l y i 


\ 


y 8 =s 2 r 2 y 2 =s 2 T 2 S 1 T l y l 


y. = S. T. 

' j 1 A 1 


i-i ^"i-i y i-i ^i-i ^i-i y i 


> 


(A10) 


y n T n ^n T n S n-i T n-i • s i T i 


J 


or v* = 


y^ = L y 
n 7 i 


where L = T S T ST 

n n-i n-i l l 


(All) 
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The boundary conditions at Station i ( Fig 20a) are 


Q l ~ k x Yt 
M i = ' K l *1 


or 


Yi ~ 5 lQi 
*i = ' CT i M i J 


( A12) 


where k, k, and 6, cr represent spring constants and influence coefficients, respectively, 
of the beam support at Station i. 

Using matrix and vector notation, both types of conditions, from equation Ai2, 
can be combined as follows: 


f 

X 

1 . . 


i 0 0 Si 


1 

i 

*1 

_ 

O 

b* 

■H 

O 


i 

M 1 


0 ■«! 1 0 


m i 

t 

tO 

t— »■ 

1 

1 

tH 

o 

o 

iH 

M 

i 






(Ai3) 


(Ai4) 


where the matrix E , is defined by equation Ai3. Only two of the four columns of S can 
be applied — the remaining two must be changed to zero columns. The possible com- 
binations of applicable columns are 


i, 2; 3, 4; i, 3; 2, 4 


(Ai5) 


What combinations should be used depends on the case in question. Obviously, 
a free end determined by 


kj = K t = 0 or f>i = °i = 

requires the combination i , 2 while a fixed end given by 
kj = K i = 00 or £>i ~ °i “ 0 

the combination 3, 4 is required. 
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If none of the four quantities k, k, 6, a, is zero or infinite any of the above com- 
binations can be selected. 


The boundary conditions of the right end can be written 


* 

Q n 

+ 

VJ 

3 

11 

o 


y n + 

6 q = 0 
n v n 



< 

- K n % “ ° 

or 

9 n - 

°n < “ 0 

> 

(A1'6) 


Using vector and matrix notation equations Ai6 change over in 


or 


1 

0 

0 

6 ~ 
n 


V 


0 

0 

1 

- a 

n 

0 


‘Pn 

= 

0 

0 

~ K n 

1 

0 


M n* 


0 

k 

n 

0 

0 

1 


V 


0 


(A17) 


: n^n = 0 


(A18) 


where E n is defined by equation A 17. 

Contrary to the case of E j , two rows only of E n can be applied while the remain- 
ing rows must be replaced by zero rows. The proper combinations are given again by 
A15. 


From the first equation All and equations A14 and A18 it follows: 

^n ^1 yi = 0 (A19) 

Since E contains two zero rows while E contains two zero columns, equation A19 
represents a linear homogeneous system of two equations. This system is solvable only 
if the determinant A equals zero. A is a polynominal in X. The roots of A = 0. 


Xj, X 2 , X 3 , . . . 


are the eigenvalues of the problem which can be found by a trial -and-error procedure. 
Insertion of Xj into equation A19 yields the solution y t , which represents the start 
vector. Now from application of equations A10 one obtains the state vectors of the 
stations and hence the i^ 1 mode shape. 
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11 


12 


(i) 


y in 

*ii 


12 


9. 

m 


The natural frequency may be calculated in accordance with equation A2. 


In the case of a free beam, equation A19 may be written: 


0 

0 

0 

0 


*11 

*12 

*13 

* 14 

0 

0 

0 

0 


*21 

*22 

CO 

CM 

*24 

0 

0 

1 

0 


*31 

*32 

*33 

*34 

0 

0 

0 

i 


*41 

*42 

00 

fc 

L 


or 


i 0 0 0 


- 

y l 


1 — 

o 

1 

0 i 0 0 


*1 


0 

0 0 0 0 
0 0 0 0 


M 1 

: Q ij 


0 

0 


* 31 Yl + * 32 <Pl = 0 
*41 yi •+ *42 <?2 = 0 


(A20) 


The frequency equation then follows as: 

* 31 * 32 _ 0 

*41 * 42 


(A21) 


Now X must be determined so that equation A2i is satisfied. With this X, equation 
A20, is solvable and yields the start vector. 




y. 

ii 




The procedure to find mode shapes and frequencies is already outlined above. 


52 



REFERENCES 


1. Bullock, T. , "Predicted Bending and Torsional Vibration of Saturn SA-D in 
Dynamic Test Tower, " IN-M-S&M-S-61-13, MSFC, Huntsville, Alabama. 

2. Liams, G. and L. L. Fontenot, "Analysis of Vibrations of Clustered Boosters, " 
ARS paper 2422-62. 

3. Storey, R. E. , "Dynamic Analysis of Clustered Boosters with Application to Titan 
HI, " AIAA Summer Meeting, Paper No. 63-208, June 1963. 

4. Loewy, R. G. , "A Matrix-Holzer Analysis for Bending Vibrations of Clustered 
Launch Vehicles." To be presented at the AIAA Structural Dynamics Symposium 
in Boston, October 1965. 1 

5. Pestel, Eduard C. and Frederick A. Leckie, "Matrix Methods in Elastomechanics, " 
McGraw-Hill, 1963. 


1. While this report was in process, the authors were informed of Loewy' s paper, which 
appears to represent a similar approach. To show the independence of this paper, the 
authors point out that their Analysis has been in use at the MSFC for two years and was 
originally submitted to the above mentioned symposium. 


NASA-Langley, 1965 M337 


53 



"The aeronautical and space activities of the United States shall be 
conducted so as to co?7tribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 

important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 

of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distri- 

bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con- 

nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



